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1 Introduction 

It is known that Normal metal - Insulator - Superconductor (NIS) tunnel junction 
in a certain range of bias voltages cools the normal metal electrode. Within a simple 
"semiconductor model" of a superconductor one can derive a cooling power of a 
single NIS junction 

P = pjjj / dE 9 -^^£M {E - eV) [f N (E - eV) - f s (E)\. (1) 

Here e is a positive absolute value of the electron charge, R is the normal state 
resistance of the NIS junction, V is the bias voltage, E is the energy of quasiparticles 
in the superconductor, A is the superconducting gap. /jv and fs are distribution 
functions in normal metal and superconductor respectively. In equilibrium these 
are Fermi functions. The cooling power (JTJ) turns out to be positive if V < A/e. 

Microrefrigerator, based on a NIS tunnel junction, has been first fabricated by 
Nahum and Martinis [TU] • They have used a single NIS tunnel junction in order to 
cool a small normal metal strip. Later Leivo et al [11] have noticed that the cooling 
power of a NIS junction (fTJ) is an even function of an applied voltage, and have 
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Figure 1: System of two NIS junctions with normal state resistances R\ and R2, Ri > ife- 
The junction number 1 cools the left normal electrode. The right normal electrode is a 
quasiparticle trap. Two possible bias methods are shown: (a) the bias voltage is applied 
between the normal metal leads, V\ 7^ 0, V2 ^ and V = Vi + V2', (b) the bias voltage is 
applied between the left normal lead and the superconductor, while the voltage drop at 
the second junction is zero, in this case V = Vi and V2 = 0. 



fabricated a refrigerator with two NIS junctions in series. By doing so they have 
achieved much better performance of the microrefrigerator. However, they have 
also observed a sharp drop of the cooling power at the base temperature below 
200 mK. This drop has been attributed to the heating of the superconducting 
electrode, which absorbs both the cooling power ([T]) and electric power VI [12] . 
Pekola et al [T3] have demonstrated that this problem can be solved if one covers the 
superconducting electrode by an additional layer of normal metal, which serves as 
quasiparticle trap and removes excited quasiparticles from the superconductor. The 
contact between the superconductor and the trap can be either direct or through 
an oxide layer. 

The aim of this contribution is to consider processes in the superconducting elec- 
trode in more detail and obtain some quantitative estimations of the effectiveness 
of quasiaprticle traps. 

The quasiparticles injected into superconductor by NIS junction create a nonequi- 
librium distribution which is characterized by the so called charge (or branch) im- 
balance. The number of electron-like excitations in not equal to the number of 
hole-like ones any more. This effect is known since 1972, when Clarke has carried 
out first experiments on charge imbalance 1 with NIS junctions. Under the exis- 
tence of charge imbalance the simple formula ((TJ) is not applicable and should be 
modified. Below we will modify it applying a method proposed by Tinkham 
Then we consider a model NISIN system with two junctions in series (Fig.l). We 
assume that the resistance of the first junction {R\) is high, while the resistance 
of the second one (R2) is low. In this case the first junction cools the left normal 
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electrode, while the second junction partially removes excited quasiparticles from 
the superconductor. The right normal metal lead models the quasiparticle trap 
[14) . For the sake of simplicity we assume that the junction number 2 is a tunnel 
junction and neglect proximity effect. 



2 Theory 

In this section we closely follow the method proposed by Tinkham [3]. 

Let us first consider a single NIS junction. For the sake of definiteness we choose 
the junction number 1. The Hamiltonian of the system can be written as follows: 



H — Hn + Hg + Ht, 



(2) 



where Hn is the Hamiltonian of normal metal, Hs is the Hamiltonian of super- 
conductor and Ht is the tunnel Hamiltonian. They are given by the following 
expressions: 

Hn = y^(e fc - eVi)c£ M Ck,cf> 

k,Ct 

Hs = ~ En + E nln, a Jn,a} , 

+ *fenK7n, Q ~ av n %,- a ]c k , a } . (3) 

Here V\ is the potential of the first normal metal lead relative to the superconductor, 
efe is the energy of electrons in the normal metal referred to the appropriate chemical 
potential, E n = + A 2 is the energy of quasiparticles in the superconductor and 
£ n are the energies of electrons in the superconductor ( has the same meaning as 
efc in the normal electrode) . The index k enumerates the states in the normal metal, 
n does the same in the superconductor, while a = ±1 is the spin index. We have 
also defined the creation c£ and annihilation c k ,a electron operators in the normal 
metal, analogous operators a^ Q and a n ,a in the superconductor, and quasiparticle 
operators 7^ a and j n ,a - The latter operators are related to the electron operators by 
means of the standard transformation a„ iQ = u fl 7n,Q — av nln -ai where the operator 
7n — a creates a quasiparticle in the state which is time reversed with respect to the 
n-th state and also carries an opposite spin. Finally, the BCS coherence factors are 



i( 



1± 



E„ 



given by the following standard expressions: |w„| 2 , \v n 

We consider the tunneling Hamiltonian as a perturbation. The current operator 
in interaction representation can be defined as the rate of tunneling out of the 
normal metal multiplied by the electron charge — e: 



Hi 



k , a 



c ta C k,OL 



(4) 
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The average value of the current can be evaluated by means of Fermi golden rule: 



(I(t))=i / dt' ([H T (t%T(t)}) 



4eRe^|i fcn | 2 J dt' 



\ Un l* e i(E n -e k+e V l){ t-t') {fStn _ fN>k) 



+ K| 2 e^" +e - eV ^- t '>(l f N . k - f s , n ) 



Evaluating time integrals we obtain 



/ = 47re ^ \tkr. 



S(E n - e k + eV)\u n \ 2 (f s ,n - /iv.fe) 



+S(E n + e k - eV)\v 2 n \(l - f N , k - f s , n ) 



(5) 



Now, as usual, we replace \t kn \ 2 by constant value |i| 2 , introduce densities of states 
in both leads, assume the distribution functions depend only on energy, and find 



h = 



1 

eR 



dt 



1 



1 + i {fs(£)-fN{E + eVi)) 



\ (l-|)CL-/s(i)-./\(-£- + M\)) 



(6) 



where 1/R = 4:ne 2 \t\ 2 N^ N§V N V S , and V N ' S are the volumes of the normal and 
superconducting leads respectively. This result has been first derived by Tinkham 
[3]. Only if the distribution function in the superconductor satisfies the particle- 
hole symmetry requirement fs(~t) — fs(t)i and the distribution function in the 
normal metal satisfies fjf(—E) = 1 — fi>f(E) Eq. ^ reduces to the standard one: 



jsy- = J_ f dE ° ^~ A2m -[f s (E) - f N (E + eVi)]. 
1 eRJ ^/E 2 - A 2 1 



(7) 



Note that there exist a difference in definition between the distribution functions 
fs(0 in @ an d fs{E) in (J7J. In the first case the distribution function depends 
on the energy of bare electrons in superconductor £, while in the second case — on 
E = ^/ 1 2 + A 2 . Moreover, negative values of E are allowed in ([7]), the distribution 
function for such values is defined by means of the relation fg(—E) = 1 — fs(E). 

The cooling power of the normal metal is obtained analogously. First we note 
that Pi = — ^2 k a £kC k a c k,cn and then repeat all the steps of the previous proce- 
dure. Thus we find 



P^ = 



1 

!Fr 



dt 



(8) 
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Again we note that only if fs(— £) = fs(0 and Jn{—E) = 1 — Jn{E) the latter 
result reduces to a simple formula (HJ. 

The same approach can be applied in order to find the rates of population and 
depopulation of quasiparticle states in superconductor. Here we already consider 
the system of two junctions. We find 

dfs,n _ d_, + . _ fs,n - fo{E n ) 

at dt To 

, \u n \ 2 fN,i{E n + eVi) + K| 2 (l - InA-Eu + eVi)) - /<?,„ 



2e 2 R 1 N^V s 

\u n \ 2 f N , 2 {E n - eV 2 ) + K| 2 (l - f N , 2 (-E n - eV 2 )) - f s , n 
2e 2 R 2 N§V s 



(9) 



where and /jv,2 are the distribution functions is the two normal metal elec- 
trodes, fo(E) = l/(exp(E/T) + 1), and r is the inelastic quasiparticle relaxation 
time in the superconductor. We use a very simple relaxation term in Eq. ([§]) be- 
cause we actually do not consider the effect of quasiparticle interactions here. The 
relaxation term is introduced in the Eq. ((3]) only in order to establish the limita- 
tions of a non-interacting model. The second term in the right hand side of Eq. ([9]) 
describes the injection of quasiparticles through the first junction, while the last 
term - the injection through the second junction. The equations similar to Eq. ^ 
are derived and discussed in Ref . [SJ [3 IS] • 

In a stationary case, /„ = 0, the distribution function of quasiparticles in the 
superconductor can be found explicitly: 

Mo - ag + _L + ^_ { j T A(E> 

1 



2e 2 i?i 
1 

2e 2 R 2 



[\u\ 2 f N ,i{E + eVi) + \v\ 2 f NA (E - eVx)] 
[\u\ 2 f N . 2 (E - eV 2 ) + \v\ 2 f N , 2 {E + eV 2 )] \ . ! ID i 



Here we have assumed that the distribution functions in the normal metal electrodes 
satisfy f^ 2 (-E) = 1 - fi 2 {E). The result KTUJ} is similar to that obtained in Ref. 
[S] , but it includes the effects of charge imbalance which become important if the 
resistances R\ and R 2 are not equal. 

The shape of the quasiparticle distribution function (I10p depends on two pa- 
rameters: the junction asymmetry parameter a — R 2 /(R\ + R 2 ) and the ratio of 
inelastic relaxation rate to the injection rate /3 = 2e 2 R\ Nq V s /tq. One can neglect 
inelastic relaxation as long as (3 -C 1 . Here we have also assumed that quasiparticle 
distribution function does not depend on space coordinates. This assumption is 
valid if the appropriate size of the superconducting electrode does not exceed the 
quasiparticle relaxation length Lj n = t/Dtq, where D is the diffusion constant. 
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Figure 2: Quasiparticle distribution function in the superconductor for different values 
of the parameters, (f) Dotted line - equilibrium distribution function /o(£) = 1/[1 + 
cxp(E/T e ff)] with effective temperature T e ///Ao = 0.3. (2) Solid line - bias voltage is 
applied between normal metals (Fig. la), eV/Ao = 1.72, T/Ao = 0.1, a = 0.05, no 
inelastic relaxation (J3 = 0). (3) Dashed line - all parameters have the same values as 
in the previous case, but the inelastic relaxation time is finite, ft = 10. (4) Short dot - 
bias voltage is applied between normal metal and superconductor (Fig. lb), eV/Ao = 1, 
T/Ao = 0.1, a = 0.05, = 0. 

Finally, the superconducting gap is determined by BCS gap equation: 



where A is the electron-phonon coupling constant. 

3 Results and discussion 

We have solved the set of equations (I6I8I10[) numerically. It has been done as 
follows. First, we have chosen the equilibrium distribution functions with the same 
temperature T in both normal metal leads. If the bias voltage is applied between 
the normal metals, (Fig. la), then the voltages across individual junctions are 
fixed by the current conservation requirement Ii(V\) = 12(^2) and V\ + V% = V. 
These equations are solved numerically. If only the junction number 1 is biased 
(Fig. lb), then we just set V\ = V, V% = 0. As long as V\ and V2 are known, we 
find the distribution function (|10p and the cooling power ([5]). Below we present 
some results. All energy scales are normalized by the superconducting gap at zero 
temperature, A = A(T = 0). 




(11) 
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Figure 3: Cooling power for different junctions asymmetry and in the absence of the 
inelastic relaxation, (a) The voltage is applied between two normal metals (Fig. la), 
(b) The voltage is applied between the normal metal and the superconductor (Fig. lb). 



Fig. 2 shows the distribution function for different sets of the parameters. In 
all cases this function has an asymmetry with respect to the chemical potential of 
Cooper pairs, this a manifestation of charge imbalance. This effect turns out to 
be relatively weak if the bias voltage is applied between the normal metals. Once 
again we note that the distribution function does not look like a Fermi function 
because of its definition (see the discussion after Eq. (J7J). 

The cooling power of the first junction ([5J as function of voltage is shown in 
Fig. 3. We have neglected inelastic processes in superconductor assuming /3 <C 
1. We observe that even a relatively small resistance of the second junction may 
significantly reduce the cooling power. It is quite natural because the cooling of the 
superconductor becomes less effective. In the opposite case /3 > 1 the cooling power 
is given by the standard formula ([1}. However, in this case the power dissipated 
in the superconductor goes first to the phonon subsystem, and then may return 
back to the cooled normal metal 12 . This effect may strongly reduce the overall 
cooling power of the device. In the absence of inelastic relaxation such processes 
are less probable, because power dissipation takes place in the normal metal with 
high heat conductivity and the heat is effectively removed from the device. That's 
why below we consider the case j3 <C 1. 

The dependence of the maximum cooling power on temperature is illustrated 
in Fig. 4. The cooling power is drastically reduced at low temperatures T/Aq < a. 
Fig. 5 shows the dependence of the maximum cooling power on the asymmetry 
a. Here we again observe the reduction of the cooling power for a > T/Aq, in 
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Figure 4: Dependence of maximum cooling power on temperature. Inelastic relaxation 
in superconductor is neglected. Solid line - maximum cooling power at a — 0.05 and bias 
voltage applied between the normal metal leads (Fig. la). Dotted line - maximum cooling 
power at a = 0.05 and bias voltage applied to the first junction (Fig lb). Dashed line - 
equilibrium result based on Eq. ©, = 0.6(Ag/e 2 i?i)(T/A ) 3/2 [Ill- 




Figure 5: Dependence of maximum cooling power on asymmetry a. j3 = 0, volt- 
age is applied to the first junction (Fig. lb). The cooling power is normalized by 
(Ag/e 2 fli)(T/A ) 3/2 . Triangles, squares and circles — numerical results for T/Ao = 
0.2, 0.1, 0.05 respectively. Solid lines — simple approximation (1121) . 
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Figure 6: The suppression of superconducting gap at different degrees of asymmetry. 
The bias voltage is applied to the normal metal electrodes (Fig. la). 



agreement with the Fig. 4. If the bias voltage is applied only to the first junction 
and the temperature is much lower than Ao, one can get an approximate analytical 
expression for the maximum cooling power: 



Pi, 



Aj /2 T 3 / 2 r 



0.76e~(^~*) a - 0.22e- 2 (^~*) Q + 0.07e- 3 (^"*) Q 



(12) 



If the voltage is applied to the normal metal electrodes, the maximum cooling power 
is even lower than (|12[) . 

Fig. 6 shows the suppression of the gap due to the heating of the superconductor. 

The model employed above requires two conditions to be satisfied: the size 
of the superconductor should be less than Lj„, and j3 should be small, /3 <C 1. 
Experimentally one finds Li n to be of the order of at least several microns (see 
e.g. |14l 113) ). Thus, the first requirement is not very restrictive, although it may 
depend on sample geometry. The parameter (3 is difficult to estimate because the 
inelastic relaxation time is unknown. Most probably this time is of the order of 
several microseconds (see e.g. [J]). If so, then the condition /3 <C 1 can also be 
easily satisfied in experiment. 

In conclusion, we have derived a generalized expression for the cooling power 
of a NIS tunnel junction taking into account charge imbalance effects. We have 
considered cooling properties of NISIN double junction structure. The interactions 
in the superconductor were neglected. Such an approximation is applicable provided 
the superconducting electrode is of submicron size. It is shown that the cooling 
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power depends strongly on the ratio of the resistances of the two junctions. We have 
shown that at low temperatures the maximum cooling power of NISIN structure is 
proportional to exp(— A R2/T(Ri + #2))- 
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